CLAIMS 

What is claimed is: 

1 . A computer program product, used in a computer system having a processor, for 
modeling an electromagnetic response of an arbitrarily shaped three-dimensional 
object to an arbitrary time-harmonic incident field by means of a well-conditioned 
boundary integral equation (BIE), the computer program product comprising: 

a, a recording medium; 

b, means, recorded on the recording medium for producing a discretized 
representation of a well-conditioned BIE to provide a well-conditioned, 
finite-dimensional linear system whose solution approximates a 
distribution of equivalent surface sources on a model caused by the 
arbitrary time-harmonic incident field; and 

c, means, recorded on the recording medium, for solving the linear system to 
determine values of the equivalent surface sources; whereby solution of 
the well-conditioned finite-dimensional linear system provides improved 
computational efficiency and solution accuracy. 

2. A computer program product as set forth in claim 1 , wherein the computer 
program product further comprises means, recorded on the recording medium for 
modeling electromagnetic scattering from the three-dimensional object due to the 
time-harmonic incident field as determined from the equivalent surface sources. 
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3. A computer program product as set forth in claim 2, wherein the computer 
program product is tailored for determining the electromagnetic scattering from a 
time-harmonic incident field in the form of a radar wavefront and where the three- 
dimensional object represents a component of a vehicle. 

4. A computer program product, as set forth in claim 1 , wherein the BIE produced is 
in a form of: 

(E excitation term)+ y(H excitation term) = {(E integral operator term)+ y{H. integral operation term 
where, 

y is a complex constant, 

J represents the unknown surface source distribution, and 
the (E excitation term) and (E integral operator term) are of a formulation 
selected from the group consisting of A and B, with 
A representing a combination of 

(E excitation term) = n(x) x E'^x)) and 

(E integral operator term) = T(k x ) T(k 2 ) , and with 
B representing a combination of 
(E excitation term) = 

tt^&X-ntt * E^OO) - C (*,)(Zn(x) • H"*(x)) and 
(E integral operator term) = 
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$kr (*,) o T T (k 2 ) + ^T\k 2 ) o T L (k 2 ) - k x k 2 T s (k,) o T s (k 2 ) ; and 

where the (H excitation term) and (H integral operation term) are of a 
formulation selected from the group consisting of A', B', and C, with 
A' representing a combination of 

(H excitation term)= Zn(x) x H mc (x) and 

(H integral operation term) = K + (k 2 ) ; 
B' representing a combination of 

(H excitation termjs (*,)(Zn(x) x H OTC (x)) and 

(H integral operation term) = K + (k^ ° K + (k 2 ) ; and 
C representing a combination of 

(H excitation term) = nx Z + (^)(Zn(x) x H mc (x)) and 

(H integral operation term) = [n x £ + (&,)]<> [n x K + (k 2 )\, 
wherein; 

T{k) ^ ikn(x) x Jdj'|G(t; x, x') J(x) + p- v[VG(fc; x, x') • J(x')]| , 
and 

s -n(x) x fa' VG(k; x, x') x J(x') , and 

s 

wherein J = Zn x H is the unknown surface current; 
£ is a complex number; 

E' nc and H'" c are incident electric and magnetic fields, respectively; 
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Z = ~\nls is the wave impedance; 

G{k; x, x') = exp(ikr)/4nr is a 3d Helmholtz kernel, 

where r = |x - x'| is a distance separating field and source points with 

harmonic time dependence e~ m assumed, and where 

T a {k)(/) = n(x)x jds'VG(k; x, x'V(x') , 

s 

T p {k)<i> s n(x) x J^n(Y ) x VG(k; x, xW) , 
s 

T\k)f = jds'VG(k; x, x'X(x') , 

s 

T T (k)f = n(x) • |<fc*VG(£; x, x')xf(x') , and 

T s (k)f = n(x)x J<&'G(A:; x, x')f(x') , and where 

s 

V and V denote differentiation with respect to the field point x and 
source point x' , respectively, and n(x)and n(x')are outward-pointing unit 
vectors normal to the target surface S at x and x' , respectively, and 

k 2 is the propagation constant in a medium external to S; 
yfcj is a complex number chosen to avoid resonances. 

A computer program product as set forth in claim 4, wherein the computer 
program product further comprises means, recorded on the recording medium for 
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modeling electromagnetic scattering from the three-dimensional object due to the 
time-harmonic incident field as determined from the equivalent surface sources. 

6. A computer program product as set forth in claim 4, wherein the BIE produced is 
in a form of: 

kT s (ik)n x E mc + iZT a (ik)n • H mc = {i [- T a (ik) o T T (k) + T\ik) o T L (k) - k 2 T s (ik) o 7* (*)]- iT (*) jj 
and wherein the terms of (ik) are developed by substituting (ik) for (k) in the terms 
of the formulation set forth in claim 4. 

7. A computer program product as set forth in claim 6, wherein the computer 
program product further comprises means, recorded on the recording medium for 
modeling electromagnetic scattering from the three-dimensional object due to the 
time-harmonic incident field as determined from the equivalent surface sources. 

8. An apparatus for modeling the electromagnetic response of an arbitrarily shaped 
three-dimensional object to an arbitrary time-harmonic incident field by means of 
a well-conditioned boundary integral equation (BIE), the apparatus comprising: 

a. a computer system having a processor, a memory connected with the 
processor, and an input/output device connected with the processor for 
receiving parameters representing an arbitrarily shaped three-dimensional 
object and an arbitrary time-harmonic incident field; 
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b. means, operating in the processor and the memory of the computer for 
utilizing the parameters representing the arbitrarily shaped three- 
dimensional object and the arbitrary time-harmonic incident field for 
producing a discretized representation of a well-conditioned BIE to 
provide a well-conditioned finite dimensional linear system whose 
solution approximates a distribution of equivalent surface sources on a 
model caused by the time-harmonic incident field; and 

c. means, operating in the processor and the memory of the computer for 
solving the linear system to determine values of the equivalent surface 
sources; whereby solution of the well-conditioned finite dimensional 
linear system provides improved computational efficiency and solution 
accuracy. 

9. An apparatus as set forth in claim 8, wherein the apparatus further comprises 
means, operating in the processor and the memory of the computer for modeling 
electromagnetic scattering from the three-dimensional object due to the time- 
harmonic incident field as determined from the equivalent surface sources. 

10. An apparatus as set forth in claim 9, wherein the apparatus is tailored for 
determining the electromagnetic scattering from a time-harmonic incident field 
represented in the form of a radar wavefront and where the three-dimensional 
object represents a component of a vehicle. 
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1 1. An apparatus as set forth in claim 8, wherein the BIE produced is in a form of: 
(E excitation term)+ y(H excitation term)= {(E integral operator term)+r(H integral operation term)}. 

where, 

y is a complex constant, 

J represents an unknown surface source distribution, and 
the (E excitation term) and (E integral operator term) are of a formulation 
selected from the group consisting of A and B, with 
A representing a combination of 

(E excitation term) = T(k x )(- n(x) x E'" c (x)) and 
(E integral operator term) = T(k } ) T(k 2 ) , and with 
B representing a combination of 
(E excitation term) = 

^r s (^)(-n(x) x E'" c (x)) - &r(k,iZn(x) ■ H"*(x)) and 

K 

(E integral operator term) = 

(*,) o T T (k 2 ) + ^r p (k 2 ) o T'iJc,) - k x k 2 T s (K) 0 T s (k 2 ) ; and 

where the (H excitation term) and (H integral operation term) are of a 
formulation selected from the group consisting of A', B', and C\ with 
A' representing a combination of 

(H excitation term) = Zn(x) x H'" c (x) and 

(H integral operation term) = K + (k 2 ); 
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B' representing a combination of 

(H excitation term)^ K\k x ){Zn{x) x IT c (x)) and 

(H integral operation term) = K + (k x ) o K + (k 2 ) ; and 
C representing a combination of 

(H excitation term) = n x K + (k x )(Zn(x) x IT C (x)) and 

(H integral operation term) = [n x K + (k^\o [n x K + (k 2 )\, 
wherein; 

T(k) s flfcn(x) x Jdy 1 |g(*; x, x')J(x) + ^ v[VG(£; x, x') • J(x')]J , 

K + (k) = ^ + K(k),md 

K{k) = -n(x) x J<fc f VG(*; x, x') x J(x') , and 

wherein J = Zn x H is the unknown surface current; 
£ is a complex number; 

E mc and are incident electric and magnetic fields, respectively; 

Z = <y////£ is a wave impedance; 

G(£; x 5 x') = exp(iAr)/47cr is a 3d Helmholtz kernel, 

where r = |x - x'| is a distance separating field and source points with 

harmonic time dependence e~ l0t assumed, and where 

T a {k)(f> s n(x) x J<fc' VG(*; x, x'^Kx 1 ) , 

^(Jfc)^ = n(x)x j* f n(x')x V*G(£; x, x>(x') ? 
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T L (k)f = fa'VG(k; x, x'X(x*), 

s 

r r (ifc)f a n(x) • jdsWG(k; x, x') x f (x') , and 

s 

T s (k)f a n(x)x Jds'G(A:; x, x')f(x') , and where 

s 

V and V denote differentiation with respect to the field point x and 
source point x' , respectively, and n(x)and n(x')are outward-pointing unit 
vectors normal to the target surface S at x and x' , respectively, and 

k 2 is the propagation constant in a medium external to S; 
is a complex number chosen to avoid resonances. 

12. An apparatus as set forth in claim 11, wherein the apparatus further comprises 
means, operating in the processor and the memory of the computer for modeling 
electromagnetic scattering from the three-dimensional object due to the time- 
harmonic incident field as determined from the equivalent surface sources. 

13. A computer program product as set forth in claim 11, wherein the BIE produced 
is in a form of: 

kT s (ik)n x E'" c + iZT a (ik)n ■ H"* = {2 ° T T (k) + T ? (ik) <> T L (k) - k 2 T s (ik) o r 5 (£)]- 

and wherein the terms of (ik) are developed by substituting (ik) for (k) in the terms 
of a formulation set forth in claim 1 1 . 
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14. An apparatus as set forth in claim 13, wherein the apparatus further comprises 
means, operating in the processor and the memory of the computer for modeling 
electromagnetic scattering from the three-dimensional object due to the time- 
harmonic incident field as determined from the equivalent surface sources. 

15. A method for modeling an electromagnetic response of an arbitrarily shaped 
three-dimensional object to an arbitrary time-harmonic incident field by means of 
a well-conditioned boundary integral equation (BIE), the method operating on a 
computer system having a processor, a memory connected with the processor, and 
an input/output device connected with the processor for receiving parameters 
representing an arbitrarily shaped three-dimensional object and an arbitrary time- 
harmonic incident field, the method comprising the steps of: 

a. utilizing the parameters representing the arbitrarily shaped three- 
dimensional object and the arbitrary time-harmonic incident field for 
producing a discretized representation of a well-conditioned BIE to 
provide a well-conditioned finite dimensional linear system whose 
solution approximates a distribution of equivalent surface sources on a 
model caused by the arbitrary time-harmonic incident field; and 

b. solving the linear system to determine values of the equivalent surface 
sources; whereby solution of the well-conditioned finite-dimensional 
linear system provides improved computational efficiency and solution 
accuracy. 
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16. A method as set forth in claim 15, wherein the method further comprises the step 
of modeling electromagnetic scattering from the three-dimensional object due to 
the time-harmonic incident field as determined from the equivalent surface 
sources. 

17. A method as set forth in claim 16, wherein the time-harmonic incident field is a 
radar wavefront and wherein the three-dimensional object represents a component 
of a vehicle, in order to determine the scattering of radar waves from a vehicle 
component due to in incident radar wave. 

18. A method as set forth in claim 1 5, wherein the BIE produced is in a form of: 

(E excitation term)+ /(H excitation term) = {(E integral operator term) + r(H integral operation term 


where, 


/is a complex, 


J represents an unknown surface source distribution, and 

the (E excitation term) and (E integral operator term) are of a formulation 


selected from the group consisting of A and B, with 


A representing a combination of 



(E integral operator term) s T{k x ) T(k 2 ) , and with 


B representing a combination of 


(E excitation term) = 
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i^ikX-nix) x E mc (x)) - ^r a (^)(zii(x) • ET(x)) and 

K 

(E integral operator term) = 

^T\k^T T (k 2 ) + ^T\k 2 )oT L (k 2 )-k l k 2 T\k,)oT s (k 2 );and 

where the (H excitation term) and (H integral operation term) are of a 
formulation selected from the group consisting of A', B', and C\ with 
A' representing a combination of 

(H excitation term) = Zn(x) x H' nc (x) and 

(H integral operation term) = K + (k 2 ) ; 
B' representing a combination of 

(H excitation term) = ^ + (^)(Zn(x) x H' nc (x)) and 

(H integral operation term)= K + (A:,) o K + (k 2 ) ; and 
C representing a combination of 

(H excitation term) = nx £ + (^)(Zn(x) x H" c (x)) and 

(H integral operation term) = [n x A: + (&,)|o [n x A" + (A: 2 )J, 
wherein; 

T(k) = ikn(x) x fa'lGik; x, x')J(x) + -L v[VG(*; x, x') • J(x')]L 


K + (k)=^+K(k),and 


K{k) = -n(x) x fa' VG(£; x, x') x J(x') , and 
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wherein J = Zn x H is an unknown surface current; 
C is a complex number; 

E mc and H " c are incident electric and magnetic fields, respectively; 

Z = -yjju/z is a wave impedance; 

G(fc; x, x') = exp(ikr)/4itr is a 3d Helmholtz kernel, 

where r = |x - x f | is a distance separating field and source points with 

harmonic time dependence e~ lG * assumed, and where 

T a (k)(f> = n(x)x J<fc'VG(*; x, x'V(x') , 

s 

T p (k)<f> = n(x) x jds'n(x')xVG(k; x, x'^(x') , 

s 

T L (k)f = faVG(k; x, x')f(x'), 

7 r (jfc)f = n(x) ■ Jdsr 1 VG(Jt; x, x') x f (x*) , and 

s 

T s (k)i = n(x)x Jdy»G(Jfc; x, x')f(x') , and where 

s 

V and V denote differentiation with respect to the field point x and 
source point x' , respectively, and n(x) and n(x') are outward-pointing unit 
vectors normal to the target surface S at x and x' ? respectively, and 
/ = V-L; 

k 2 is the propagation constant in a medium external to S; 
k x is a complex number chosen to avoid resonances. 
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19. A method as set forth in claim 18, wherein the method further comprises the step 
of modeling electromagnetic scattering from the three-dimensional object due to 
the time-harmonic incident field as determined from the equivalent surface 
sources. 

20. A method as set forth in claim 18, wherein the BIE produced is in the form of: 
kT s (ik)n x E mc + iZT a (ik)n ■ H' wc = jf [- T a (ik) o T T (k) + T*(ik) o T L (k) - k 2 T s (ik) o T s (k)]- K + (k)}j 

and wherein the terms of (ik) are developed by substituting (ik) for (k) in the terms 
of the formulation set forth in claim 18. 

21 . A method as set forth in claim 20, wherein the method further comprises the step 
of modeling electromagnetic scattering from the three-dimensional object due to 
the time-harmonic incident field as determined from the equivalent surface 
sources. 
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